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Abstract 

We study reversible quantum cellular automata with the restriction that these 
are also Clifford operations. This means that tensor products of Pauli operators 
(or discrete Weyl operators) are mapped to tensor products of Pauli operators. 
Therefore Clifford quantum cellular automata are induced by symplectic cellu- 
lar automata in phase space. We characterize these symplectic cellular- automata 
and find that all possible local rules must be, up to some global shift, reflection 
invariant with respect to the origin. In the one dimensional case we also find 
that every uniquely determined and translationally invariant stabilizer state can 
be prepared from a product state by a single Clifford cellular automaton timestep, 
thereby characterizing these class of stabilizer states, and we show that all ID 
Clifford quantum cellular automata are generated by a few elementary opera- 
tions. We also show that the con^espondence between translationally invariant 
stabilizer states and translationally invariant Clifford operations holds for peri- 
odic boundary conditions. 
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1 Introduction 

Classical cellular automata have become a standard modeling tool for complex 
phenomena. With their discrete time step and their intrinsically high degree of paral- 
lelization they are ideally suited for models of diverse phenomena as coffee percola- 
tion, highway traffic and oil extraction from porous media. As an abstract computa- 
tional model cellular automata can simulate Turing machines, and even explicit simple 
automata such as Conway's life game have been shown to support universal compu- 
tation yj. On the quantum side, the interest in cellular automata stems from their 
implementation in optical lattices and arrays of optical microtraps. However, the the- 
ory of quantum cellular automata (QCAs) is still in its early stage. Since each cell may 
influence several others, the dynamics is subject to a "no-cloning" constraint, leading 
to a non-trivial interplay between the conditions of locality and unitarity. 

It is therefore helpful to have some class of QCAs, which can be analyzed in great 
detail, and which can serve as a testing ground for general ideas about QCAs. This pa- 
per is concerned with such an analysis, namely of the special class of Clifford quantum 
cellular automata (CQCAs), in which the elementary time step is given by a "Clifford 
gate", meaning that it takes tensor products of Pauli matrices to tensor products of 
other Pauli matrices. In the theory of gate model computation, and for the one-way 
quantum computation model, a detailed analysis of what can be done with Clifford op- 
erations alone turned out to be very useful, even though - as the downside of allowing 
an efficient classical description - such gates alone do not allow universal quantum 
computation. By analogy it is therefore clear that CQCAs do not comprise the full 
complexity of QCAs. What one can hope to get, however, is an interesting class of 
cellular automata, and some tools for understanding this class in great detail. 

A similar analysis has been done with Gaussian quantum cellular automata e.g. 
the QCA describes a chain of harmonic oscillators with nearest neighbor couplings. 
For all these QCAs the Hilbert space of one elementary cell is infinite dimensional, 
and the QCA maps phase space translations, also referred to as Weyl operators, to 
phase space translations. In our approach we use elementary cells with a finite number 
of levels, which corresponds to replacing the continuous phase space by a discrete 
space. 

1.1. Definition of Clifford quantum cellular automata. — By definition, a cellular 
automaton is a lattice system, which consists of many subsystems (called "cells") la- 
beled by a point lattice in space. For simplicity, we will always take the lattice as Z'^, 
the integer cubic lattice in s space dimension^. The cell systems in the classical case 
may have states like "occupied" and "empty". In the quantum case, they will be p- 
state quantum systems, for some finite p. In either case the group of lattice translations 
("shifts") is a symmetry of the system. 

The dynamics will be given by a discrete global time step, or global "transition 
rule" assumed to have the following three properties: 

'see however Section]?] where we discuss periodic boundary conditions, and hence toroidal lattices 
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• translation invariance: the time step commutes with the lattice translation sym- 
metries. 

• reversibility: there is an inverse rule. For a finite quantum system this would 
mean unitary dynamics. For the infinite lattice system this will be stated alge- 
braically below. 

• locality, or "finite propagation speed": the state of each cell after one step can 
be computed from the state in a fixed finite region around the cell. 

These assumption define the class of reversible QCAs O. The locality and reversibil- 
ity conditions are best phrased in the Heisenberg picture: if A denotes some observable 
of the system, its expectation after one time step starting from the initial state p will be 
{T{A))p for a suitable observable T{A), where by {A)p we denote the expectation of 
A in the state p. The transformation A ^ T{A) is what we will call the global rule of 
the automaton. Then reversibility (together with complete positivity, which is required 
of any dynamical map) implies that T is a homomorphism of the observable algebra of 
the whole system: T is linear, T{AB) = T{A)T{B), and T{A*) = T{A)*. Locality 
means that an observable A^ localized at some lattice point x ^If (i.e., an observable 
of the cell at x) will be mapped to an observable localized in the region x + M . That 
is, T(Aa,) will be in the tensor product of the cell algebras belonging to the sites x + n 
with n G M . By translation invariance this set A/", called the neighborhood scheme of 
the automaton is independent of x. 

The global transition rule T is a map on an infinite dimensional space, and hence 
not readily specified explicitly. However, by using the basic properties of QCAs one 
can see that it suffices to know just a few local data, associated with the region A/", in 
order to reconstruct T uniquely. Suppose we know T(/l^) for every observable A^ in 
some basic cell x. Then by translation invariance we know the analogous transforma- 
tion for any cell. Moreover, since every local observable can expanded in products of 
one-cell observables, and because T is a homomorphism, we can compute T for any 
local observable. So the restriction T^. of T to the observables of a single cell x can be 
called the local rule of the automaton. We can also decide by a finite set of equations, 
whether a proposed local rule actually belongs to a well-defined global rule: clearly 
Tx must be a homomorphism. The only further condition one has to check is that the 
images of observables A^ and By localized in the cells indicated also commute, i.e., 
Tx{A^Ty{By) = Ty{By)Tx{Ax), whcncvcr x ^ y. This is necessary, because A,j. and 
By commute, and a moment's reflection shows that this is also sufficient for uniquely 
reconstructing the images of arbitrary observables under T. The commutation con- 
ditions on the local rule are trivially satisfied when x and y are sufficiently far apart, 
when X + H and y + M are disjoint. Hence only finitely many conditions need to be 
checked. 

For special classes, the job of specifying a QCA via its local rule can be reduced 
still further, which is where the Clifford condition comes in. Let us assume now that 
we have a qubit system, so the local cell dimension isp = 2. For each local cell we thus 
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have a basis for the observables, consisting of the identity and the three Pauli matrices, 
which we denote by X, Y and Z. By etc. we denote the corresponding Pauli matrix 
belonging to the cell x. Finite tensor products of Pauli matrices belonging to different 
sites, perhaps with a sign ±1, ±i will be referred to as Pauli products. These form a 
group, called the Pauli group. Then a Clifford quantum cellular automaton (CQCA for 
short) is defined by the condition 

• Clifford condition: If A is any multiple of a Pauli product, so is T{A). 

Clearly, this is equivalent to the property that one-cell Pauli operators are taken to Pauli 
products, which simplifies the local rule. Moreover, it suffices to specify T(X^) and 
T{Zx) for some x, because we can compute T{Y.j) = TiiX^Zx) = iT{Xx)T{Zx) 
via the homomorphism property. Hence a CQCA is defined in terms of just two Pauli 
products. 

Example 1.1. For the one-dimensional lattice (s = 1), consider the relations 

TiyXx) = Zx 

T{Zx) = Zx-i® Xx ®Zx+i ^ ^ 

Let us verify that all requirements for a local rule are satisfied. To begin with each 
of the expressions on the right hand side, as a product of Pauli matrices, is hermi- 
tian with square one. These are all the required conditions related to just a single 
line, and are satisfied for any Pauli product with a sign ±1. Next we have to ver- 
ify the anti-commutation relation arising from applying a homomorphism T to the 
anti-commutation relation XZ + ZX = Indeed, T(X^)T(Z^) + T(Z^)r(X^) = 
-Zx-i ® (Z^Xx + X^Zx) ® Zx+1 = 0. Hence the definition T(F^) = iT{Xx)T{Zx) 
again produces a hermitian operator with square 1, and the local rule is a homomor- 
phism into the algebra on the sites x + J\f with J\f = { — 1, 0, 1}. Finally, we have to 
check the commutation rules for the images of observables from neighboring sites. For 
example, we have [T{Xx),T{Z,j.+i)] = —[Z^, Z^ ® X^+i ® Z,j.+2] = 0, and similarly 
[T{Z.j.),T{Zx+2)] = 0. Perhaps the only non-trivial relation to check is 

[T(Z,), T(Z,+i)] = ® X, ® Zx ® Xx+i ® Zx+2] = 0, 

which holds because the factors on sites x and x + 1 both anri-commute. 

In principle, we would also have to check the existence of an inverse for the au- 
tomaton, which is actually given by T(Xx) = X^^i ® Z^. ® X^+i and T{Zx) = —X^, 
but as was shown in |[2l, this already follows from the homomorphism property. 

It is clear from this example that the search for CQCAs is now a combinatorial 
problem. We can first look for self commuting Pauli products, i.e., Pauli products, 
which commute with all translates of itself. Only these can appear on the right hand 
side of local rules. One can then check, for any pair X', Z' of such products, whether 
they anti-commute, while all proper translates of X' commute with Z'. In fact, we 
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began our investigation by running this simple search program. We found, for exam- 
ple, that while there is a rich variety of self-commuting Pauli products only reflection 
symmetric products could appear in a local rule. This will indeed be shown in full 
generality below. 

1.2. Translationally invariant stabilizer states. — Commuting sets of Pauli products 
also play a central role in the problem of determining so-called stabilizer states: these 
are pure states, which can be characterized by eigenvalue equations for Pauli products 
or, equivalently, by the condition that certain Pauli products have expectation ±1. It 
is easy to check that Pauli products which simultaneously have sharp expectations ±1 
must commute. Now for the infinite lattice systems it is natural to ask which Pauli 
products A have the property that there is a unique pure state p of the infinite system, 
which has expectation 1 for A and all its translates. 

As the simplest example, let us take A = Z^, so we ask for states with (Zx) = 
1 for all X E Z'^. Clearly, this defines the "all spins up" state, which is an infi- 
nite product state. A slightly more complex example uses the stabilizer operators 
A = Zx^i Xx Zx+i, which singles out the one-dimensional cluster state, whose 
higher dimensional analogs are used as the entanglement resource for universal one- 
way quantum computing jSl. 

Showing that these eigenvalue equations define a unique state of the infinite lattice 
is now very easy, by using the cellular automaton ([T]): Since this automaton maps Z^ to 
the required stabilizer operator, all existence and uniqueness problems for such a state 
are mapped to the corresponding trivial questions for the stabilizer operator Z^. In 
other words, self-commuting Pauli products of the form A = T{Zx) for some CQCA 
T characterize a unique translation invariant cluster state. We will show later that (at 
least in one dimension) the converse is also true, so that there is a very close connection 
between stabilizer states and Clifford cellular automata. 

1.3. Our methods and techniques. — The definition of CQCAs given above applies 
only to qubit systems. However, all our results are also valid for higher dimensional 
cells, particular cells of prime dimension p. The role of the Pauli operators X and Z is 
then taken by the cyclic shift on C^, and the multiplication by a phase, i.e. 



where all ket labels q are taken modulo p. Products of these operators are called Weyl 
operators, and the appropriate definition of CQCAs requires that T{Xr^) and T^Z^) 
are both tensor products of Weyl operators. The necessary preliminaries on the Pauli 
group and Clifford operations in this extended setting, and the background concerning 
infinite lattice systems are provided in Subsection l2.1[ 

In order to utilize the translation symmetry one would like to use Fourier transform 
techniques. However, in the discrete structures an integral with complex phases makes 
no sense. It turns out, however, that a "generating function" technique does nearly as 
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well. The analogue of the Fourier transform is then a Laurent-polynomial in an inde- 
terminate variable, i.e., a polynomial with coefficients in the field Zp = Zp with both 
positive and negative powers. The salient facts about this structure will be provided in 
Subsection l2.3[ 

The description in terms of Laurent polynomials can also be adapted to lattices 
with periodic boundary conditions. This will be described in SectionlH 

1.4. Outline and summary of results. — In order to discuss general Clifford quan- 
tum cellular automata, that is, for arbitrary lattice and single cell dimension, we in- 
troduce in Section [2] the necessary mathematical tools. We first review the concept of 
discrete Weyl systems and (infinite) tensor products of them, thereby characterizing 
the underlying "phase space". We show that Clifford QCAs can be completely char- 
acterized in terms of classical symplectic cellular automata. We also introduce our 
Fourier transform techniques and study the structure of isotropic subspaces, because 
these play an essential role for the characterization of symplectic cellular automata and 
translationally invariant stabilizer states. 

In Section [3] we will state our main results. We show that symplectic cellular 
automata can be identified with two-by-two matrices, which have Laurent-polynomials 
as matrix elements. We will find that the determinant of this matrix must be one and 
that the polynomials must be reflection invariant. In the one-dimensional case we state 
that every translationally invariant stabilizer state can be prepared out of a product state 
by a single CQCA step. Furthermore, we also specify the generators of all ID QCAs. 

Finally, we show in Section |4] that the close connection between translationally 
invariant stabilizer states and CQCAs also holds in the case of periodic boundary con- 
ditions even in every lattice dimension. 



2 Mathematical tools 

We introduce some mathematical tools, which we will use to study Clifford QCAs. 
We start with a short repetition of finite Weyl systems, which generalize the Pauli 
operators to systems with prime number dimensions. These Weyl operators can be 
described by phase space vectors and Clifford operations are induced by symplectic 
transformations on the phase space. Since we are looking for translationally invariant 
operations, we also introduce some kind of Fourier transform. 

2.1. Weyl algebras. — Each single cell in a QCA is given by a finite dimensional 
quantum system, so the observables on a single system can be described by matrices 
from the algebra Aip{C). A possible basis for this algebra is given by Weyl operators 
w(r, k) = X''' Z'^, whereby X and Z are given by the generalized Pauli operators from 
equation Q. These operators fulfill the Weyl relations 

w(ri + ra, h + A^a) = e-''^^'w{n, h)w{r2, fcs) , (3) 
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where e = exp(27ri/p) is the p^^ root of unity. ^From this equation the commutation 
relation 

w(ri, h)w{r2, k2) = e'''^'-''^'w{r2, k2)w{n, h) (4) 
immediately follows. Obviously we get for p = 2 the standard Pauli operators from 

X = w(l,0), r = zw(l,l), Z = w(0,l), (5) 

and the Weyl operators are generalizations of the Pauli operators to higher dimensional 
spaces. The indices r and k are integers modulo p, so they are elements of the finite 
field ¥ = Zp. In infinite dimensional systems Weyl operators describe phase space 
translations and therefore we call the space a discrete phase space. 

Building a tensor product of Weyl operators means that we must assign a phase 
space vector ^(x) = {^+{x),^_{x)) E to each lattice point x G Z^, so is a 
mapping from Z'^ into F^ and we denote for the tensor product 

w(0 = (g) w(e(x)) . (6) 

This infinite tensor product is well defined, if there are only finitely many of the Weyl 
operators different from 1 = w(0). For a mapping ^ : Z'^ ^ F^ we have that only 
finitely many x with ^(x) ^ are allowed, so the support of ^ is finite. The set of such 
functions describes the global system and is identified with the global phasespace Eg. 
We denote the finitely supported functions from to F by Cf(Z*) and we have Eg = 
Cf{Z^Y. The corresponding Weyl operators generate an algebra and, by restricting 
the support of the functions to some finite subset A G Z"*, we get a finite dimensional 
algebra ^(A) = (S)a;gA ■Mp{C), also called the local algebra of A. By taking the union 
of these algebras over all finite subsets of Z* and taking the closure (in operator norm) 
we get a quasilocal C*-algebra A (H, which is used in the general theory of quantum 
cellular automata 121 . 

The local structure is accompanied by the symmetry group of lattice translations. 
For each lattice translation x G Z'^ an automorphism is defined by 

«x.w(0 = w(r,.0 . (7) 

where is the translation of phase space vectors. Given a phase space vector ^, the 
translated vector is {t^^ ){y) = ^{y — x). So the automorphism a.j. shifts the position of 
each tensor factor by x. It follows directly from (|7]) that the homomorphism property 
ax+y = axUy holds. Furthermore, the automorphism a.^ maps the local algebra .4.(A) 
onto ^(A + x). 

The Weyl relations of a single system completely determine the relations of the 
global system which are given by 
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where we have introduced the bilinear form 77) := Ylix&^ ^+i^)V-{^)- The ad- 
joint of a Weyl operator is given by 

w(0* = e"^«'«)w(-0 (9) 

which is due to the unitarity of the Weyl operators. 

Since commutation relations are essential for validating possible local rules of 
quantum cellular automata, the commutation relations of Weyl operators are most im- 
portant for us. We get 

w(r/)w(0 = e''(^''')w(Ow(r/) , (10) 

whereby cr{^,ri) := /3{^,ri) — (3{ri,^) is the canonical symplectic form on Eg. This 
means that two Weyl operators w(^) and w(?7) are commuting if and only if a{^, rj) = 
(and for p = 2 they anti-commute if (t(^, 77) = 1). In particular, an abelian algebra of 
Weyl operators is given by a subspace of Eg on which the symplectic form vanishes. 
Such a subspace is called isotropic and a maximally abelian algebra corresponds to a 
maximally isotropic subspace. 

2.2. Clifford quantum cellular automata. — As already mentioned a Clifford quan- 
tum cellular automaton is a QCA which maps Weyl operators to multiples of Weyl 
operators, which are in our case tensor products of single cell Weyl operators, so we 
have the relation (the "Clifford condition") 

nw(o) = ^(Ow(te) (11) 

with a mapping t on the phase space Eg and some phase valued function ^9 : — > 
W(l) = {z E C\\z\ = 1}. Since T is an automorphism we find with equation (flOl) that 
w(t^)w(t?7) = e°'^^'^^w{tri)w{t^) holds, so we have a{t^,tr]) = cr{^,r]) or in other 
words t is a symplectic transformation. 

For reversible operations the Clifford condition is in general equivalent to the Weyl 
covariance (for general theory on covariant channels we refer to [7] and for the special 
case of Weyl covariance to |[8l|9l) of the quantum channel: 

Proposition 2.1. An automorphism T on the Weyl algebra A fulfills the Clifford con- 
dition 477]) if and only if the Weyl covariance 

T{w{r])Aw{r])) = w(tr/)r(A)w(tr/)* Vr/ G Eg (12) 

holds for all operators A E A and some symplectic transformation t. 

Proof. Because the Weyl operators form a basis of A we just have to insert w(,^) for 
some ,^ G in the covariance condition, which yields the equation e'^^^'^'^T{w(^)) = 
w(t?7)T(w(,^))w(t?7)*. If T is a Clifford automorphism we have already seen that 
t is a symplectic transformation and obviously fulfills this equation. In the inverse 
direction we get that T{w{^)) must be a multiple of w(t(^), because the relation must 
hold for all r] E Eg and the symplectic form is non degenerate (note that the support 
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of the phase space vectors is finite and that t maps therefore finitely supported vectors 
to finitely supported vectors, so the commutation relations can be checked in a finite 
dimensional space). □ 

Since a QCA is a translationally invariant automorphism on the quasilocal alge- 
bra, it suffices that the Clifford condition holds for the local rule, e.g. the QCA re- 
stricted to operators which are localized in a single cell. Furthermore, because of the 
Weyl relations on a single cell, we only need to specify the image of the Weyl op- 
erators w(l, 0) and w(0, 1). To some extend we are free in the choice of the phases 
■!9(l,0),'i9(0,l) G W(l), since these phases do not interfere with the commutation re- 
lations for the local rule. The only condition is that some power of a Weyl operator is 
always equal to 11 (we will specify this below), and so these phases must be some roots 
of unity. The two phases 0) and i3{0, 1) completely determine the function 'd. 

Of course t and ^9 must be translationally invariant, because T is translationally 
invariant. Using the homomorphism property of the QCA and equation ([8]) we get 
+ r?)w(t(^ + 7])) = ?9({)^?(r/)e'3(«''')-'3(*«'*'?)w(t^ + trj), so - because the Weyl 
operators form a basis - the transformation t must be linear and the phase function 
must fulfill 

^{^ + r/) = ^9(^)|?(r;)e/^«''?)-/5(t5.t^) , (13) 

which enables us to calculate the phase ?9(^) for each ^ E Eg, if the local rule and 
therefore t and the phases '(9(1, 0) and ^?(0, 1) are given. In total we get the following 
theorem: 

Theorem 2.2. If T is a Clifford quantum cellular automaton ( equation ([77]) ) on the 
Weyl algebra A, then t is a translationally invariant linear symplectic transformation 
("symplectic cellular automaton" ) and the phase function fulfills equation ( 17?] ). 

This means that we are able to study Clifford QCAs - up to some phase function 
- in terms of a classical cellular automaton on the phase space Eg. It is well known 
that Clifford operations allow an efficient classical description, which in the case of 
QCAs turned out to be the group of classical symplectic cellular automata. In the 
rest of the paper we will study the structure of this kind of cellular automata, thereby 
characterizing the structure of CQCAs. 

We would like to give a closed expression for the phase function 19, but this has to 
be done in dependence of the cell dimension. First we consider the case p ^ 2. Then 
all Weyl operators fulfill w(^)p = 1 and because of T{t) = 1 the phase 'd{^) must 
be a p*^ root of unity. So we can write ^?(^) = e*^*^^^ with a function Lp : Eg ¥, 
which then has to fulfill + rj) = (p{^) + ip{ri) + (3{^, rj) — /3(t^, trf). This equation 
determines the function 93 (^) up to some linear functional A(^), which is given by 
the choice of the phases 0) and (y9(0, 1). The bilinear form I3{^, rj) — P{t^, trj) is 
symmetric, because t is a symplectic transformation. If ]? 7^ 2 we may divide by 2 and 
the general solution is ip{0 = - tO) + KO- 

The case of qubits (p = 2) is slightly more complicated because the Weyl operators 
fulfill w(0^ = {-if^^'^H. So the phase function must fulfill ^9(0 = i^^^^ with if : 
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Ss ^ Z4. We replace the form : x — > Z2 by the bilinear form (3 : Eg x Eg ^ 
Z4, which is formally given by /3 = 2/3, so the values of P are even elements of Z4 
and the Weyl relation becomes w(,^ + v) = i'^^^''^^w{^)w{r]). This means that if fulfills 
ipi^ + V)= ^(0 + V^(^) + 7(^, V) with the form 7(e, r]) = t]) - l3{t^, tr/). This 
form is symmetric, so in the decomposition 7(,^, 77) = J2i j lij^iVj we have 7^^ = 'jji 
and all these elements are even. We can find 7j with 'ja = 2'ji, but this choice is not 
unique in Z4 and corresponds exactly to the freedom in the choice of the phases 0) 
and 'd{0, 1). The solution for (/? is then given by v?(0 = Y.i<j lij^i^j + Y.i lS (note 
that e {0, 1} and so = holds). 

2.3. Algebraic Fourier transform. — We would like to use Fourier transform tech- 
niques for the study of the structural properties of symplectic CA, because of trans- 
lational invariance, and because we know that this is very helpful for symplectic CA 
with continuous single cell phase space [3]. So we have to apply a Fourier transform 
to the functions Cf{7j^ ). But the values of these functions are in the finite field F and 
multiplying such a value with a complex number does not really match. It turns out 
that a slight modification of the usual Fourier transform does as well. For a function 
/ e Cf(Z^) we define 

with = ■ ■ ■ Ug" . Now the transformed function / is a polynomial or, more pre- 
cisely, a Laurent-polynomial in the variables ui, . . . ,Us with coefficients in F, which 
will be denoted by Vg = F[mi, ■ ■ ■ , u^^, ■ ■ ■ , m^^]. Note that we have indeed poly- 
nomials, because the functions in Cf(Z'^) are finitely supported. Equation (fT4l) iden- 
tifies functions of Cf{Z^) with polynomials Vg and this identification is unique, so 
Cf(Z^) and Vg are isomorphic. The usual Fourier transform would require u„ = e*^". 
We do not further specify the domain of the variables, and this approach can be seen 
as "generating function approach" or "algebraic Fourier transform". 

The convolution /t^t/i = J2x f{~^)'Txh is a natural produc|§of functions in Cf(Z'*). 
The invertible elements with respect to this operation are the functions which are sup- 
ported on a single lattice point, e.g. / = c6x (Sx is the Kronecker-delta) with c G F 
and X G Z", and the unit element is 60. The nice fact about Fourier transform is that 
the convolution turns into a usual product which is also true for our algebraic version: 

f7h = fh f,heCf{Z'). (15) 

Note that the invertible polynomials are monomial^ e.g. they are of the form u^. 
Of course the unit element is the constant 1 = 5o. Another important operation is 
the reflection operation (or involution) f{x) := f{—x) for / e Cf(Z,*). Obviously 
the reflection preserves the convolution, e.g. f ^ h = f ^ h, and for the transformed 

function we have f{u) = f{u~^). 

-With the convolution the set Cf(Z'') becomes a "commutative division ring". 
^This will be different when we go to periodic boundary conditions. 
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The phase space Eg consists of two-dimensional tuples of functions from Cp{Z'^) 
and all operations can be defined component-wise0, so we get that the phase space is 
isomorphic to Eg = Vg. We would like to study the structure of symplectic CA in 
this polynomial space. The transformation of an operation t : ^ is defined 
according to t,^ = t,^, so t is a mapping from Vg to Vg. We introduce the symplectic 
form a-.V^x Vl Vg by 

a(e,r/) =e7^--e^r/+ , (16) 

which can be written as (t(^, 77) = det(^, 77), whereby 77) denotes the 2 x 2-matrix 

{i.v)=[f_ :;:) (IV) 

with polynomial entries. The symplectic form a is the best fitting symplectic form for 
symplectic CA, because it combines both the basic symplectic form a as well as the 
translation invariance: 

Proposition 2.3. A linear operation t on the phase space "Eg is a symplectic cellular 
automaton, if and only if, the transformed operation t leaves the symplectic form cr 
invariant. 

Proof. For this proof we introduce the form a{^,r]) = a{^,T(.)r]) for ^,r] E Eg. A 
straightforward computation shows that a{^, 77) = ^+ i^r]^ — ^_ t^tt/^ holds. This means 

we have 5'(^, 7/) = (t(^, 77), so a is the Fourier transform of a and the invariance of a 
under some operation t is equivalent to the invariance of a under t. 

Now suppose t is a symplectic CA. Then we have for all x Elf that ^(t,^, tr/) (x) = 
cr(t{, TjXt]) = a{t^, tTxT]) = a{^, TxT]) = a{^, v){^) holds, because t is translationally 
invariant and preserves a, so a is invariant under t. 

If t leaves a invariant, this holds also for a = cr(-,-)(0), and because of this 
(j(t^, r^.t7/) = a{^,Txr]) = (T(t^, tr^-T/) holds for all x E If and all ^,7/ E Eg and 
so t must commute with the translations Tn^. □ 

So we can characterize symplectic CA in "momentum space" by studying the linear 
transformations on Vl which leave the symplectic form invariant. In the subsequent 
we will mainly work in the polynomial space Vg. Therefore we will just identify the 
phase space Eg with and we will omit the symbol " for the Fourier transform of 
transformations. 

2.4. Isotropic subspaces. — As we have already seen in the introduction, commu- 
tation relations are important for the verification of local rules of reversible QCAs, 
because a QCA is a homomorphism and preserves the algebraic structure. Especially 
the images of and must be "self-commuting", meaning that \r{X,^,T{Xy)\ = 



With the component-wise convolution the phase space is a two-dimensional CF(Z'')-module. 
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= [T(Zx),T{Zy)] holds for all x,y E If . So the operators T{X,j) generate a transla- 
tionally invariant abelian algebra. For Weyl operators translationally invariant abelian 
algebras correspond exactly to isotropic subspaces of with respect to the symplectic 
form a and these subspaces can be easy connected to translationally invariant stabilizer 
states. Therefore it is important for us to study the structure of these subspaces. 

A Ps-subspace|f|x C is called isotropic, if for all 77 G X the symplectic form 
ct[^, 77) = vanishes. An isotropic "P^-subspace X is called maximally isotropic, if the 
relation ?]) = for all ,^ G X implies that 77 G X holds. 

For us the form of the generators of isotropic, in particular maximally isotropic, 
Ps-subspaces is important, because this is a substantial step for the characterization 
of local rules of CQCAs and translationally invariant stabilizer states. The following 
lemma shows that a generator ^ of a singly generated maximally isotropic subspace 
is reflection invariant and that the components ^+ and are coprime. We will call a 
polynomial f E Vg (or a tuple of those) reflection invariant for some half integer lattice 
point a E ^Z^*, if / = n^"/ holds. The greatest common divisor of two polynomials 
f,h E Vs will be denoted by gcd(/, h). Note that the greatest common divisor is 
defined only up to invertible elements. We will simply write gcd(/, h) = 1, if / and h 
are coprime. 

Lemma 2.4. 

1. If the subspace Vs^ C V'^ is maximally isotropic, we have gcd(,^+, ^J) = 1. 

2. If the subspace Vsi C is maximally isotropic, ^ is reflection invariant to 
some point a E 

3. Every reflection invariant polynomial generates an isotropic Vs-subspace. 

Proof. Ad 1. Suppose Vsi is a maximally isotropic P^-subspace and gcd(^_|., = / 
is not invertible. So we can write ^ = frj with gcd(?74., rfj) = 1, but 77 ^ Vs^ since / is 
not invertible. But we have that cr(^, 77) = fcr{^,^) = holds, which is a contradiction 
to Vs^ being maximally isotropic. 

Ad 2. Suppose that Vg^ is maximally isotropic. By 1 we have gcd(^+,.^_) = 1. 
Since a(^, ^) = - ^_^+ = 0, it follows that ^ G Vs^. So we have ^ = with 
some polynomial f E Vs- But for the reflected phase space vector ^ we also have that 
gcd(^+, ^_) = 1, so / must be invertible and therefore a monomial f{u) = for 
some a E ^Z^. 

Ad 3. Suppose that ^ = u^'^^ is reflection invariant. Then ^) = (t(m^"^, ^) = 
— i-ijr) = holds, and ^ generates an isotropic Ps-subspace. □ 

Example 2.5. Both ^1 = (1 + u)i^^ and ^2 = reflection invariant. The 

corresponding Weyl operators w(^i) = Zq® Zi and w(^2) = -Z^-i ® Xq® Zi are the 
same reading from the left and from the right ("palindromes"). Both phase space vec- 
tors generate isotropic subspaces. The subspace generated by ^2 is indeed maximally 



^More precisely one should say submodule, but we will use the more convenient word subspace. 
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isotropic and the components ^2,+ and ^2 - are coprime, whereas the subspace gener- 
ated by ^1 is not maximally isotropic because 1 + m is a nontrivial common divisor. 
This is also clear in terms of operators, because all operators commute with w(,^i), 
but cannot be obtained by products of translates of w(^i). 

In particular, the greatest common divisor comes into play. We will be able to 
state more results in the one-dimensional case (s = 1), due to the fact that the ring 
of polynomials V := V\ = ¥[u,u^^] is euclidear|^. Especially this means that the 
euclidean algorithm can be applied for finding the greatest common divisor of two 
polynomials, which is also used for the factorization of wavelet transformations ifTSl . 

Lemma 2.6 (Extended euclidean algorithm for Laurent polynomials). Let ^ E be 

a phase space vector. Then there exist /o, /i eV such that 

/oe+ + /i^ = gcd(e+,e-) (18) 

holds. 

Proof. We define the degree of a Laurent polynomial / = Yl^ti. f^'^^ ^^^U) •= 
L+ — L_ when /i_ and /^^ are nonzero. Suppose deg(^+) < deg(^_) and let tq = 
,^4. and ri = We make a division with remainder and get a polynomial go with 
deg(go) = deg(ro) — deg(ri) and a polynomial r2 with deg(r2) < deg(ri) such that 

To = qon + r2. (19) 

With this decomposition we get gcd(ro,ri) = gcd(ri,r2). We repeat this division 
recursively until the remainder vanishes: 

Ti = qiVi+i + ri+2 (20) 

Tn+l = Qn+lTn- (21) 

Then we have r„ = gcd(r„, r„+i) = gcd(ro, ri). We rewrite the recursion to get the 
form of equation (fTSl ): 

' r^-l A 1 \f r,_2 
Ti ) I 1 -qi-2 ) \ n-i 



So we get 



; r„...ro( 



with 

r 



1 

1 -q^ 

and since all entries in the matrices are polynomials we get polynomials /o and /i such 
that 

Tn = foro + flTi 

holds. □ 



^In more abstract words is a principal ideal ring, which means that every ideal in V is generated 
by a single element. For this general algebraic theory we refer to ifTOl . 
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3 Main results 

3.1. Characterization of Clifford quantum cellular automata. — We have seen in 
Proposition |23] that symplectic cellular automata are nothing else but linear functions 
on the phase space Eg = that preserve the "Ps-symplectic form a. Such a map on 
Vg can be represented by a two-by-two matrix with entries in the polynomial ring Vg. 
The first column is given by ti = t (p) ("the local rule for X") and the second column 
by t2 = t d") ("the local rule for Z"). The commutation relations of the local rule then 
end up in the following conditions on the column vectors: 

Corollary 3.1. A two-by-two matrix t with entries in Vg is a symplectic cellular au- 
tomaton, if and only if, the column vectors of t = {ti,t2) fulfill a(ti,ti) = = 
(j(t2,t2) and a{ti,t2) = 1. 

Remark 3.2. The column vectors ti 2 of a symplectic cellular automaton generate 
maximally isotropic "P^-subspaces "Psti 2, since these are the images of the basis vec- 
tors (q) and (°) under the invertible symplectic transformation t. Because the basis 
vectors generate by construction maximally isotropic subspaces, this must then also be 
true for the images ti^2- 

In the next subsection, we shall see that the classification of one-dimensional 
symplectic cellular automata is easier to handle. A useful observation is that a s- 
dimensional symplectic cellular automaton t G ^^2(^5) induces for each direction 
k = 1, . . . , s a one-dimensional cellular automaton. To see this, we introduce for each 
direction k = 1, . . . , s a surjective ring homomorphism which maps the polynomial 
ring Vg of s-variables ui, - ■ ■ ,Ug, u^^, ■ • • , Ug^ onto the ring V of one variable u. The 
ring homomorphism assigns to a polynomial f E Vg the polynomial 

r,/(n):= . . . , x^) n^'" (22) 

(x-i,...,a;'')eZ'' 

which only depends on the variables u, u~^. The ring homomorphism evaluates the 
polynomial f E Vg at ui = 1, for / =^ k, whereas Uk = u is the remaining free variable. 

For a symplectic cellular automaton t E M.2{Vg) the conditions (^(ti 2,ti 2) = 
and (T(ti,t2) = 1 are identities of polynomials. The matrix r^t E M.2{V) is build 
by applying the ring homomorphism to each matrix element individually. Obvi- 
ously, the identities rfc(j(ti 2, ti 2) = o'(rfc'ti,2) rfcti^2) = as well as rfe(j(ti,t2) = 
(T(rfcti, rfct2) = 1 follow. As a consequence we get: 

Corollary 3.3. Let t E M.2{V's) be a s-dimensional symplectic cellular automaton. 
Then for each direction k = 1, . . . , s, the two-by-two matrix r^t E 7^2 (P) <^ one- 
dimensional symplectic cellular automaton. 

Now it is easy to show that symplectic cellular automata are reflection invariant and 
that the determinant is a monomial. It is slightly more involved that we have reflection 
invariance with respect to a lattice point and not with respect to an half integer lattice 
point. 
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Theorem 3.4. A Vs-Unear map t G J^2{Vs) is a symplectic cellular automaton, if 
and only if, the following holds: 

1. The matrix t is a reflection invariant with respect to some lattice point a G Z**. 

2. The Vs-valued determinant oft is det(t) = m^". 

Proof. If t is a symplectic cellular automaton, then the column vectors ti 2 generate 
maximally isotropic subspaces. By Lemma [24l it follows that ti, respectively t2, is re- 
flection invariant to some half integer lattice point a, respectively h. Since t preserves 
the symplectic form a we obtain 1 = (T(ti, = (j{u^'^ti, u^^i2) = u'^^^~"'^a{ti, t2) = 
^2(fe-a) therefore a = b for an half-integer lattice point a G ^Z**. As a conse- 
quence, t is reflection invariant for a G ^Z'^. Now, 1 = a(ti,t2) = o-(u^°ti,t2) = 
det(ti,t2). 

Vice versa, let t be a matrix, which is invariant with respect to the reflection at a 
and whose determinant is det(t) = u'^"-. Then the column vectors ti 2 are reflection 
invariant, which implies (by Lemma that (5"(ti^2, ti,2) = holds. The determinant 
of t is which implies a(ti, t2) = 1. Thus t preserves the symplectic form a. 

By Corollary |3.3[ we obtain a one-dimensional symplectic cellular automaton r^t 
for each lattice direction k = 1, . . . , s. We have already shown that the column vectors 
ti 2 are reflection invariant for a = (a\ . . . , a^), which implies that for each direction 
k the column vectors r^ti 2 are reflection invariant for a*^. We also have that r^ti 2 
generate maximally isotropic P-subspaces, since these define valid cellular automaton 
rules. 

Suppose now, that / G P is reflection invariant for 6 G |Z in the half-integer 
lattice. Then we can translate / by an even translation 2y G 2Z, such that c = 2(b + y) 
is either or 1. If / is of even length, then c = 1 follows. The polynomial u^f is 
reflection invariant for 1/2 and can be expanded as 

Now, for each n E N, the polynomial + is a multiple of m + 1. Thus / 
is also a mu ltiple0 of M y{u + 1). From this we conclude that, if b is not an integer, 
then a reflection invariant ^ G is a multiple of u~y{u + 1) and does not generate 
a maximally isotropic P-subspace, since n + 1 is a nontrivial common divisor of 
and which is a contradiction. So must be an integer lattice point, that is, a = 
{a\...,a') eZ'. □ 

So each symplectic cellular automaton t is reflection invariant for the reflection 
at some lattice point a G Z"*. Therefore, the symplectic cellular automaton M~"t is 
reflection invariant with respect to the origin a; = 0. In the subsequent, we call all 
symplectic cellular automata, which are reflection invariant with respect to the origin. 



^Note that the coefficients are from the finite field Z. 
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to be "centered" and it is sufficient to classify only those. The polynomials in Vs 
which are reflection invariant with respect to the origin form a subring TZg C Vs and 
will be simply called reflection invariant (for s = 1 we will again omit the index). 
From Theorem 13 .41 we obtain a handy characterization of centered symplectic cellular 
automata: 

Corollary 3.5. The group of centered symplectic cellular automata is given by the 
group SL(2, TZg) of two-by-two matrices t with entries in the subring TZg of reflection 
invariant polynomials and TZg-valued determinant det{t) = 1. 

Example 3.6. The symplectic transformation corresponding to the "cluster state QCA" 
(eq. 11.11) is given by 

■ 1 

1 u + u'^ 



(24) 



Obviously all entries are reflection invariant with respect to the origin and the determi- 
nant is equal to one (modulo 2). 

Remark 3.7. A nice aspect of Corollary 13.51 is that the centered symplectic cellular 
automata can be obtained by the following strategy: Choose two arbitrary reflection 
invariant f,hETZs and find all possible factorizations of the polynomial fh — l = f'h' 
into a product of two reflection invariant /', h' E IZg. The corresponding symplectic 
cellular automaton is then given by 

I, GSL(2,7^,). (25) 

Even if the task of factorizing the polynomial fh — \ h quite cumbersome, there is 
always a "trivial" solution, namely, h' = 1 and f = fh — 1. The matrix 

( ) ^SL(2,7^,) (26) 

describes the corresponding symplectic cellular automaton. 

Remark 3.8. Another remarkable fact is that, due to Cramer's rule, the inverse of a 
centered symplectic CA t is simply given by 

-1 / t22 — tl2 



t-^ = . (27) 

\ —^21 111 / 

Similarly we have that for a symplectic CA t containing a translation by a positions, 
e.g det(t) = u"-, the inverse contains a translation by —a positions. 
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3.2. ID CQCAs and translationally invariant stabilizer states. — In this subsec- 
tion we are investigating one-dimensional symplectic cellular automata. As already 
mentioned, we can achieve more results in this case, because we can apply the eu- 
clidean algorithm (Lemma [Z6l) . We will use the euclidean algorithm to show that for 
every reflection invariant ^ E V'^ with gcd(.^+, = 1, there exists at least one cor- 
responding reflection invariant rj such that d-{r], = 1 holds and is therefore a valid 
column of a symplectic cellular automaton matrix. We will use this fact to show that 
every uniquely determined and translationally invariant stabilizer state can be prepared 
from a product state by applying one timestep of a Clifford QCA. 

Stabilizer states are studied extensively in the last years ( IfTSl and [[T6l are just 
examples, which are useful as introductory texts). The basic concept is to fix an abelian 
group of operators (usually a subgroup of the Pauli group), also called stabilizer group, 
and to define a stabilizer state as common eigenvector of all these operators. In our case 
we are looking for translationally invariant states, so the stabilizer group is generated 
by all translates of one single Weyl operator w(^) for some phase space vector ^ G V^. 
The state should fulfill uj{axw{^)) = uj{w{tx^)) = 1 for all x E Z. The stabilizer 
formalism is often studied for finitely many qudits. In that case it is known that the 
stabilizer state is uniquely determined, if the number of generating operators is large 
enough (see e.g. |fT6l for a quantitative statement). In our situation we have infinitely 
many qudits, so we cannot apply this result. But it turns out that the operators w(rj,^) 
must generate a maximal abelian algebra, or equivalently, the subspace must be 
maximally isotropic. 

Theorem 3.9. For a phase space vector E, eV^ the following is equivalent: 

1. There exists a uniquely determined state uj with uj{w{txC)) = ^far all x G Z. 

2. VC, is a maximally isotropic V-subspace. 

3. There is a Clifford QCA T with w(0 = T(w(0, 1)). 

4. C, is a reflection invariant and gcd(^+, ^_) = 1. 

Proof. [2l =^ m Because is a maximally isotropic subspace we conclude from 
Lemma [24] that ^ is reflection invariant with gcd(^+, ^„) = 1. 

m ^>|3l We have to find t] E with d-{r],^) = 1 and a(r], rj) = 0. With Lemma 
|2.6l we find a solution f±EVof the equation /+^+ + /-C- = gcd(^+, ^_) = 1 and 
rj' := (/_, — /+) is a solution of d-{r]',^) = 1. Yet we do not know, whether 7]' is 
reflection invariant, or equivalently, whether cr{'r]', r]') = holds. But if r]' is a solution 
of (t(?7', ,^) = 1 then the same is true for rj = rj' + f^. Thus we have to solve the 
condition = a(r/, r]) = a{ri', r]') + a{f^, 7]') + a{ri', f^) = a{r]', rj') - J +f - The 
polynomial h := a{'r]', r]') is anti-symmetric with respect to the reflection ^ ^ ^ and it 
can be expanded as /i = X]n>o ~ u~^). By choosing / = X]n>o ^nU^ we find 
that 7] = r]' + is indeed reflection invariant. The matrix 77) G SL(2, TZ) is then a 
symplectic cellular automaton and induces a Clifford QCA with the desired property. 
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|3l Consider a state uj with the desired property. Then this state is equal to u)o 

T, whereby is a state with a;(a;a;w(0, 1)) = 1 for all x G Z, so the stabilizer group of 
this state is given by all translates of w(0, 1). This means is a translationally invariant 
product state, which is determined by the equation cj(w(0, 1)) = 1 and corresponds 
to the one dimensional projector onto the eigenspace of w(0, 1) with eigenvalue 1. 
Therefore this state is uniquely determined and uj = o T is the unique state with 

[U =^ |2l Suppose is an isotropic P-subspace but not maximally isotropic. 
By Lemma IATI we know that there exists a phase space vector ^ with ^ . 
So we have ,^ = with / not invertible and Lemma [Z4l tells us that ^ is reflection 
invariant. With help of the euclidean algorithm we find a QCA T and a corresponding 
symplectic transformation! withT(w(0, 1)) = w(t(0, 1)) = w(,^') (just as step two of 
this proof). Now consider a product state with (/?(w(0, u^)) = exp{^ax) depending 
on the a^. We transform this state with and the expectation values of the operators 
Tzw(^) should be all equal to 1: 

1 = (/p^-i(r,w(0)=^(r,.w(t-iO) 

= ^(r,.w(t-i(/0)) = ^(r.w(/^ (M))) 

= ^(w(5^/_,4+x(0,l))) =n/-'i^(w('5/c+.(0,l))) 

k k 




k 



So we have to solve the equations J2k f-k^k+x = to get appropriate and therefore 
states with the desired property. Since / is not invertible the support of / is not a 
one-elementary set. Let / = {—L_, . . . , —Lj^} be the minimal interval such that 
supp(/) C /. We can choose arbitrary ai_, . . . ,aL+-i to compute a^^ from the 
equation Ylik f-kO-k- Recursively all can be calculated from the other equations but 
the solution will depend from the initial choice of the aL_ , . . . , a^^.i. This means that 
there exists more than one state (p of the above form, such that 1 = (/^^--i (r^.w(,^)) is 
fulfilled. So the uniqueness of the state in[TJ forces to be maximally isotropic. □ 

So we have shown that every translationally invariant and uniquely determined 
stabilizer state in a one-dimensional lattice can be prepared out of a product state by 
a single timestep of a Clifford QCA. Unfortunately we cannot generalize this result to 
higher lattice dimensions, because Lemma [Z6] is only valid for univariate polynomials. 
The euclidean algorithm for computing the greatest common divisor can be general- 
ized to multivariate polynomials |fT9ll , but the extended version (equation (fTSl) ) does 
not hold. 

Example 3.10. We consider again the phase space vectors = (1 + u)(^^) and 

^2 = (u+u-i) (see Example 12.51) . As already mentioned, the phase space vector 
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is reflection invariant for 1/2 and generates an isotropic P-subspace, but none of the 
statements of Theorem 13 .91 holds: The expectation value of w(tx^i) = Z^+i is 
equal to one both in the "all spins up" and in the "all spins down" state, so there is 
no uniquely determined stabilizer state. As we have seen in 12.51 the subspace V^i is 
not maximally isotropic. The reflection invariance does not hold for an integer lattice 
point, so ,^1 is not a valid column of a symplectic CA, and 1 + n is a common divisor 
of ,^1,+ and which is not invertible. 

In contrast ,^2 fulfills all four conditions. The uniquely determined stabilizer state 
is given by the one-dimensional cluster state and a possible CQCA is given by exam- 
ple O 

3.3. Factorization of ID Clifford QCAs. — We have seen that the set of centered 
CQCAs form a group and that this group is given by 2 x 2-matrices with determinant 
one and reflection invariant polynomials as matrix elements. In the one-dimensional 
case the group structure can be more clarified, since we are able to give a complete set 
of generators, which can be regarded as elementary operations. 

A simple example of a 2 x 2-matrix in SL(2, 7^^) is for some reflection invariant 
polynomial f eTZ given by 

g(/) ■■=(] ? V (28) 



which we will call "shear transformation". In particular, g(/i + = g(/i)g(/2) 
holds for all /i,/2 G 7?I§. The symmetric polynomials Wn = + n G N, 
and Wq = 1 form a basis of the subring TZ. Thus every shear transformation can 
be decomposed into a finite product of elementary shear transformations g{cWn) with 
n G N U {0} and c e F. 

The local rule of the corresponding QCA Gn with G'„w(^) = w(g(w„)^) is for 
n > 1 given by 

GniZpi) = Zq 

For n = we have the single cell operation ("local shear transformation") 

Go(Xo) = w(l,l) 

GoiZo) = Zo ' ^^""^ 
which correspond for p = 2 to applying the phase gate ( ^ . ) to all single cells. 



Another single cell operation is the "local Fourier transformation", which is in 
phase space given by the matrix 



'^This means, the map / 1-^ g(/) is a group homomorphism from the additive group TZ into the 
group of centered symplectic cellular automata SL(2, TZ). 
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with some constant 7^ c G F (for c = 1 we will write f := fi). For p = 2 we have 
c = 1 and the corresponding QCA switches the operators X and Z in each single cell 
and is therefore given by applying the Hadamard matrix. 

Note that all symplectic single cell transformations can be obtained by a product of 
local shear and local Fourier transformations, which is a generalization to higher cell 
dimensions of the fact that local Clifford operations are generated by Hadamard and 
phase gate. 

The symplectic transformations g(/),fc are elementary symplectic cellular au- 
tomata in the sense of the following theorem. The proof, which is technically slightly 
more involved, is given in the Appendix |A.2[ 

Theorem 3.11. The group of centered symplectic cellular automata SL(2, IZ) is gen- 
erated by the set {g(w;„)|n G N U {0}} U {fc|c G F}. 

Remark 3.12. A more concrete formulation of the statement of Theorem 13. 11 1 is that 
every one-dimensional centered symplectic cellular automaton t is a finite product of 
shear transformations and local fourier transforms of the following form: 

t = g(/.)fc.---g(/2)fc.g(/l)fc, (32) 

with reflection invariant polynomials fi, . . . , fr E TZ and constants ci, . . . , G F. 

Example 3.13. Let us consider in the qubit case (p = 2) the symplectic cellular au- 
tomaton 

Wi 1 
I + W2 Wi 

and, since W2 = wi"^ for p = 2, we have det(t) = 1. The corresponding CQCA is 
given by 

T{Zo) = Xo ®Zi ■ ^^^^ 

The basic idea for deriving a decomposition like (l32l) is to reduce the support of the 
first column of t by applying a shear transformation from the right. We get 

tg(u,) ={° I)- 05) 

This matrix is obviously equal to g{wi)i and we have 

t = g{wi)fcg{wi) , (36) 

which is indeed a decomposition in accordance with (|32|) . 

Remark 3.14. For p = 2 all the generators g(/) and f are their own inverses, so the 
time evolution of these operations alternates between the identity and a single timestep 
of the automaton. Especially these QCAs show no propagation, because the neigh- 
borhood of the iterated automaton does not increase with the number of timesteps. A 
nontrivial time evolution only occurs, if the symplectic cellular automaton is composed 
of at least two different generators. 



(33) 
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4 Periodic boundary conditions 

In this chapter we are looking for translationally invariant Clifford operations with 
periodic boundary conditions on an s-dimensional lattice. These boundary conditions 
are given by an s-dimensional torus T|^, which is determined by s independent lat- 
tice vectors N = {Ni, . . . , Ns) (see figure [1]), and all lattice points which differ by 
these vectors are identified. The number of (not identified) lattice points is given 
by I T^l := | det(A''i, . . . , Ns)\. We denote here by Vs^n the ring of polynomials 
/ = XlxGT^ f{x)u^ such that the n-variables fulfill the periodic boundary conditions 
= = . . . = = 1. This guaranties that the product of two polynomials 
from Vs,N is again an element from Vs,n- But algebraically there are large differences 
between Vg and Vs,n' 'Ps,n is not a division ring, because there are zero divisors and 
there are in general other invertible elements than = u^^U2^ ■ ■ ■ The reflection 
/ is again given by replacing u by u^^ or in other words we substitute by u^~^. 
The symplectic form a is then of the same form as in the infinite lattice case. 




Figure 1: A 2-dimensional torus defined by A^i = (1, 3) and N2 = (5, 1). 

Now we have to say, how a Clifford QCA (or a symplectic cellular automaton) is 
defined on these systems. In the general theory of QCAs the neighborhood of a 
QCA with periodic boundary conditions is not allowed to be too large in comparison 
with the torus. This guaranties that the QCA can be extended to the whole lattice. 
Since this case is covered by restricting the existing Clifford QCAs to periodic bound- 
ary conditions, we drop all locality conditions, and we take the same structure as in 
Corollary 13. H as definition of a symplectic cellular automaton: 

Definition 4.1. A 2 x 2 matrix t = (ti, with entries in Vs,n is a symplectic cellular 
automaton if the column vectors fulfill a{ti, ti) = = a(t2, and (T(ti, t2) = 1. 

With this definition it is possible to state an analogous version of Theorem 13 .91 also 
for periodic boundary conditions. But the proof is quite different from the infinite 
lattice case. 

Theorem 4.2. For a phase space vector eV^ the following is equivalent: 

1. There exists a uniquely determined state uj with uj{w{t.j.^)) = Ifor all x G T|^. 

2. Vs,Ni is a maximally isotropic Vs,N-subspace. 
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3. There is a symplectic cellular automaton t with ^ = ty^. 

Proof. [B <^=^|2l A stabilizer state on M qudits is uniquely determined, if and only if, 
the minimal number of generators of the stabilizer group equals M |fT6l [T71 . Here we 
have the | | generators w{t,j,E) ■ These are independent, if and only if, they generate a 
maximally abelian algebra, or equivalently, if Vs,nC, is a maximally isotropic subspace. 

121 ^^[3l Since we have a finite dimensional space, there exists a symplectic basis, 
and any basis of a maximally isotropic subspace can be extended to a symplectic basis 
ll20l . For this construction we turn to the original phase space (by inverse "Fourier 
transform") and define by := r^^ basis vectors of a subspace. Since we know 
by [2l that this space is isotropic, these vectors fulfill o-{^x,^y) = and therefore 
= ^) = (t(^, Then there exists a dual vector fj with a{fi, = ^xo and we 
define rjy = Tyfj. We get that (j{'qy,ix) = (^iv^^x-y) = ^xy holds and so we have 
^{v^O = 1- We have to verify that we can choose r], such that a{r], r]) = holds. But 
if ?7 is a solution to a(?7, = 1 the same is true for f]' = r] + for any / G Vs^n and 
we can find similar to the case of Theorem 13.91 an appropriate / with a{r]', r]') = 0. 

|3l =^ [21 Suppose that t is a symplectic cellular automaton with ^ = t d*) . Then 
t induces a homomorphism between the maximally isotropic subspace Vs,n (i) = © 
Vs,N and Vs,N^ with t/Cj*) = t(°) = Since t is invertible and preserves the 
symplectic form a, it follows that any maximally isotropic subspace is mapped onto a 
maximally isotropic subspace, which implies that Vs^n^ is maximally isotropic. □ 



Example 4.3. As an example of a translationally invariant stabilizer state, we consider 
translationally invariant graph states ifTTl . The graph is encoded by its adjacency ma- 
trix r = (r(x, y))a;,a;eT^, and the isotropic subspace that determines the graph state is 
given by the phase space vectors 



r/ 



G = Cv{T%f (37) 



with / G Cf(T^). Translation invariance of the graph state implies that the matrix 
elements r(x, y) depend only on the difference x—y, so there is a function 7 G Cf(T^) 
such that r(x, y) = 7(x — y) holds. Thus after Fourier transform the phase space 
vector, which generates the maximally isotropic subspace, is given by ,^ = (]^) and we 
can choose a suitable symplectic cellular automaton t with ^ = t (°) by 



1 7 
1 



(38) 



Figure [^represents a translationally invariant graph state on the ID torus Tg = Zg. 
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Figure 2: A translationally invariant graph state on the ID torus Tg = Zg. 



The adjacency matrix T is given by 



v 



1 

1 1 
1 1 1 
1 1 



1 



(39) 



and obviously only depends on the difference x — y of the variables x,y E Zg. Apply- 
ing the Fourier transform, yields the polynomial ^ = + + = + + u^. 
The symplectic cellular automaton, which creates the graph state as explained above 
is given by the matrix 



t 



1 u + + 







Note that t is reflection invariant, since 



1 

is a reflection fix-point. 



(40) 




We are going to present another example to show that the phase space vectors do 
not have to be reflection invariant, because not all invertible elements are monomials. 
But the invertibility of a fixed polynomial depends on the size of the torus and so a fixed 
phase space vector ^ may define a translationally invariant stabilizer state for some N, 
but it is possible that there exists A^', such that Vi^N'i is not maximally isotropic, and 
therefore ^ does not characterize a unique stabilizer state for A^'. 



Example 4.4. We consider for p = 2 the phase space vector 



^ = (1 + u + v!^) 



u 



u 



(41) 



on a one dimensional torus of variable size. Note that = holds for all A^, so ^ 
generates an isotropic subspace. The corresponding tensor product of Pauli operators 
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is given by 



(42) 



and is obviously not reflection invariant. 

Let us first have a look at the case = 7. It is easy to show that 1 + u + is not 
invertible. We define f = {/(I + u + u^) = (""'^+") and have that = 0, but 

^ ^ Vij^. So Vij^ is not maximally isotropic, and there is no unique stabilizer state. 

For = 6 the inverse of 1 + n + n'^ is given hy u + u"^ + u^, so ^ and ^ generate the 
same subspace, which is actually maximally isotropic. So ^ is indeed a valid column 
of a symplectic automaton, but starting from the "all spins up" state both CQCAs 
corresponding to 



5 Conclusions 

We have analyzed the structure of Clifford quantum cellular automata that act on 
a s-dimensional lattice of p-level systems. The results which can be achieved depend 
on the dimension of the lattice and whether we put periodic boundary conditions or 
working with the infinite lattice. 

We have characterized the group of CQCAs in terms of symplectic cellular au- 
tomata on a suitable phase space. With the help of Fourier transform, this phase 
space can be identified with two-dimensional vectors of Laurent-polynomials, and 
symplectic cellular automata can be described by two-by-two matrices with Laurent- 
polynomial entries. We have shown that these entries must be reflection invariant and 
that up to some global shift the determinant of the matrix must be one, so the group of 
CQCAs is isomorphic to the special linear group of two-by-two matrices with reflec- 
tion invariant polynomials as matrix elements. 

We have proven that there is a correspondence between ID CQCAs and ID trans- 
lationally invariant stabilizer states. For a fixed translationally invariant pure stabilizer 
state u^^, which is in particular a product state, every other translationally invariant 
pure stabilizer state if can be created by applying an appropriate CQCA to the 
chosen product state: if = u;®^ o T^. 

Pure stabilizer states can be also characterized by maximally isotropic subspaces. 
We have characterized the phase space vectors, which generate maximally isotropic 
subspaces, namely their components must be coprime and reflection invariant. 

In the one-dimensional case we have also more clarified the group structure of 
CQCAs. As we have shown, each one-dimensional CQCA can be decomposed into 




prepare the same stabilizer state. 
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a product of elementary shear automata and local Fourier transforms, so the group of 
CQCAs is generated by this set of operations. 

For periodic boundary conditions the techniques from infinitely extended lattices 
can be applied to a certain extend. According to the discussion of translationally invari- 
ant stabilizer states on the ID lattice, we have proven that there is an analogous corre- 
spondence between CQCAs and translationally invariant stabilizer states with periodic 
boundary conditions even in any lattice dimension. 
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A Proofs and technicalities 

A.l. Ad Theorem 13.91 — For the proof of Theorem 13 .91 we need that a singly gen- 
erated isotropic subspace can always be embedded into a singly generated maximally 
isotropic subspace: 

Lemma A.l. Let ^ & and Vsi be an isotropic, but not maximally isotropic 
Vs-subspace. Then there exists a phase space vector rj E V"^ such that Vsf] ^ Vs^ is 
maximally isotropic. 

Proof. Vs^ is isotropic if and only if the equation = a{$,, ^) = ■C+'C- ^ ^-^+ holds. 
We make a distinction of cases for this equation: 

i. ^+ = (analogously ^_ = 0): Then Vs^ = {0} © Vs^- and ^_ is not invertible 
since this would force Vs^ = {0} © Vs to be maximally isotropic. So we can set 

77 =(0,1). 

ii. ^+ = f^^ (analogously = f^+) with / reflection invariant: Then we have 
Vs^ = Vsi-{f, 1). We set 77 = (/, 1) and get that Vs'T] is a maximally isotropic 
subspace since = a{ri, A) = /A_ — implies A = A_?7 G VsTj. 

iii. ^+ 7^ 7^ ^_ and ^± ^ f^^: Then ^ = with / invertible, so { is reflection 
invariant for some n E ^Z'^. Because Vs^ is not maximally isotropic we can find 
T] ^ Vs^ withO = d-{^,ri) = n"({+?7_— ,^„?7+). Since and,^_ are nonvanishing 
this implies ^ = grj for some g E Vs- We can choose gcd(?7+, rfJ) = 1 and VsT] 
to be maximally isotropic. 

□ 
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A.2. Ad Theorem 13.111 — For a polynomial f E V the coefficient of the monomial 
is {f)x- Recall that "degree"of a polynomial in f E V is, defined by deg(/) : = 
max{x\{f)x 7^ 0} — min{a;|(/)a; ^ 0} and that the support is defined by supp(/) : = 
{x\{f).^0}. 

Lemma A.2. Let rj) be a symplectic cellular automaton which is invariant under 
the reflection at the origin: ^ = ^ and rj = f]. If the degrees of column vectors 
fulfil deg(.^) > deg(?7) then there exists a shear transformation g(/), with reflection 
invariant f eV, such that the symplectic transformation 

(e',V) = (e,^)g(/)fi (43) 

satisfles deg{^,T]) > deg{^', rj') and deg{^') > deg{r]'). 

Proof. Since ^ and t] are reflection invariant, the degree is an even integer and we 
introduce x := deg(,^)/2, y := deg(?7)/2, as well as rii := x — y > 0. We conclude 
from the identity cr{^,r]) = 1 that 

(a(e,r/)),+, = {^+)^Ar]~)y " {^-)~Ar]+)y = (44) 

is valid. This implies that 

(e). = (e)-. = -/i(r/), = -/i(r/)-, (45) 

for some /i E ¥. Now + + u~"-^)r]) = which implies that 

deg(e + + n""^)r/) < deg(0 • (46) 

Now we observe 

(6,r/i) := (e,r/)g(/iK^+«-'^^)) 

from which we conclude that deg(^i,?7i) < deg(^,?7). If deg(,^i) > deg(?7i) we can 
find a shear transformation g{f2{u^^ + u~^^)) such that 

(6,^2) = (ei,r;i)g(/2K^+n-"^)) (48) 

fulfills deg(^25 V2) < deg(^i, ?7i). We can proceed this reduction until the /th step with 
21 = |supp(.^) \ supp(?7)|. The resulting symplectic cellular automaton 

{^i,Vi) = (6-1, ^^-i)g(/K«"' (49) 
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then satisfies deg(^/) < deg(?7/). If deg(^i) = deg{rji) then, there is an appropriate 
constant fi+i E F such that 

(r, V') = (e, ^)g(/)fi = i-Vi, + fi+iVi) (50) 
holds with deg(.^') > deg(?7'). Here / is the reflection invariant polynomial 

i+i 

If deg(^z) < deg(?7z), then the shear transformation g(/;+i) is not applied and we get 

(e',V) = (e,^)g(/)fi = (-^ao (52) 

with the polynomial / = Yl]=i fji'^"'^ + ^ 

Proof of Theorem \3.11\ Let (^07^70) be a symplectic cellular automaton which is in- 
variant under the reflection at the origin. Then, by Lemma lAT2l there exists a symplec- 
tic cellular automaton (^1, 771) and a shear transformation g(/i) such that 

(eo,^o) = (6,^i)g(/i)fi (53) 

and deg(,^i) > deg(?7i). Thus we can iterate this reduction process until {^k, Vk) is a 
constant symplectic transformation (corresponding to a QCA with single cell neighbor- 
hood), which can be decomposed into a product of local shear transformations g(wo) 
and local Fourier transforms fc- with q G F. This yields the following decomposition 

of (^o,'7o): 

(eo, VO) = g(/r)fc. ■ ■ ■ g(/2)fc.g(/l)fci . (54) 

□ 
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